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Abstract 

Approximate analytical solutions of the Boltzmann equation for particles 
that are either extremely relativistic or non-relativistic when they decouple 
from the thermal bath are well established. However, no analytical formula 
for the relic density of particles that are semi-relativistic at decoupling is yet 
known. We propose a new ansatz for the thermal average of the annihilation 
cross sections for such particles, and find a semi-analytical treatment for cal- 
culating their relic densities. As examples, we consider Majorana- and Dirac- 
type neutrinos. We show that such semi-relativistic relics cannot be good cold 
Dark Matter candidates. However, late decays of meta-stable semi-relativistic 
relics might have released a large amount of entropy, thereby diluting the den- 
sity of other, unwanted relics. 
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1 Introduction 



The accurate determination of cosmological parameters by up-to-date observations, 
most notably by the Wilkinson Microwave Anisotropy Probe (WMAP) [1], increases 
the importance of quantitative predictions. In particular, the estimate of the cosmo- 
logical relic abundances of particle species is essential, for the history of the universe 
depends on these quantities. One of the most important examples is the cosmological 
abundance of dark matter [21 [3] , whose mass density is found to be [Ij 



f2 is the mass density in units of the critical density. 

The abundance of some particle species is determined by solving the Boltzmann 
equation, which describes the change of the particle number caused by particle reac- 
tions as well as by the expansion of the universe [2]- However, there is no analytical 
general solution of this nonlinear differential equation, and therefore one needs to 
solve the equation numerically in many cases. In early studies, approximate analyt- 
ical formulae have been found for the relativistic [H [2] and non-relativistic [51 El El 
regimes. 

Stable or long-lived weakly interacting massive particles (WIMPs) with weak- 
scale masses are examples of cold relic particles, which decouple from thermal equi- 
librium when they are non-relativistic. In standard cosmology, decoupling of WIMPs 
occurs in the radiation-dominated (RD) era after inflation, and analytic approximate 
formulae for the WIMP relic abundance have been derived [5l [6l [7j. In the opposite 
limit, where decoupling of particles from the thermal background occurs when they 
are relativistic, the relic abundance is approximated by its equilibrium value at the 
decoupling temperature, and not sensitive to details of its freeze-out [2J . The result- 
ing relic density can therefore easily be computed analytically. On the other hand, 
no analytical treatment to calculate the relic abundance in the intermediate regime 
is known yet. 

In this paper, we revisit the relic density of particles x that decouple from 
the thermal bath when they were semi-relativistic, i.e. at freeze-out temperature 
Tp ~ m^. Assuming that the Maxwell-Boltzmann distribution can be used for all 
participating particles, we introduce an expression for the thermal average of the x 
annihilation cross section which smoothly interpolates between the extremely rela- 
tivistic and non-relativistic regimes. It is shown that our new ansatz is capable of 




1 



reproducing the exact thermally-averaged annihilation cross section with accuracy 
of a few percent. Given this approximated cross section, we can define the freeze-out 
temperature by comparing the annihilation rate to the expansion rate. The assump- 
tion that the comoving x density remains constant after freeze-out turns out to be 
a good approximation for the relic abundance of semi-relativistic particles. 

We also discuss the roles such semi-relativistic particles could play in realis- 
tic cosmological scenarios. It should be emphasized that the abundance of semi- 
relativistically decoupled relics tends to be large because it is only very mildly Boltz- 
mann suppressed. We point out that scenarios where semi-relativistically decoupling 
particles form the Dark Matter have problems with structure formation, Big Bang 
Nucleosynthesis (BBN) and/or laboratory measurements. Nevertheless, such semi- 
relativistic relics can be useful for diluting the density of other, unwanted relics by 
late-time out-of-equilibrium decay [SiiSj. As an example, we investigate a scenario of 
decaying sterile neutrino that is assumed to depart from thermal equilibrium when 
it is semi-relativistic, in sharp contrast to non-thermal sterile neutrino scenarios 
|10j . Thermal equilibrium is attained by introducing some higher-dimensional oper- 
ator. It is illustrated that an enormous amount of entropy can be produced without 
spoiling the successful BBN prediction of the light element abundances. 

This paper is organized as follows: In Sec. 2 we begin by reviewing briefly the 
method to calculate relic densities for relativistic and non-relativistic particles. In 
Sec. 3 we explain the new formalism which is applicable for all freeze-out tempera- 
tures in case of S- and P-wave cross sections. The way to calculate the freeze-out 
temperature is also shown. In Sec. 4, the possibility for semi-relativistic particles to 
have the observed dark matter relic density of ^dm^^ ~ 0.1 is considered. Then, we 
discuss the amount of entropy produced by the decay of unstable semi-relativistic 
species that decay in less than a second. Finally, Sec. 6 is devoted to summary and 
conclusions. Some properties of modified Bessel functions are described in Appendix 
A. In Appendix B we argue that the use of the Maxwell-Boltzmann distribution in 
the definition of the thermally averaged cross section only leads to a small mistake 
in the final relic density. 
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2 Relic abundances in the non— relativistic and rel- 
ativistic limits 



In this Section we briefly review the standard analytical approximations for evalu- 
ating the relic abundance of hypothetical particles x [21 El [7]. These are applicable 
to particles that were either non-relativistic or extremely relativistic at freeze-out. 

The number density is obtained by solving the corresponding Boltzmann equa- 
tion [2]. For the moment, we assume that single x production and x decay are for- 
bidden by some symmetry or adequately suppressed. The Boltzmann equation takes 
a simple form if one further assumes that the quantum statistics factors describing 
the Rose enhancement or Fermi suppression of all final states can be neglected; in 
Appendix B it is shown that this is essentially equivalent to assuming that the distri- 
bution functions of all relevant particles are proportional to the Maxwell-Boltzmann 
distribution. The Boltzmann equation for can then be written as [2] 



where ?T-x>eq is the x equilibrium number density, {av) is the thermal average of 
the annihilation cross section multiplied by the relative velocity between the two 
annihilating x particles, and H is the Hubble expansion rate of the universe. The 
second term on the left-hand side describes the dilution caused by the expansion of 
the universe; the first (second) term on the right-hand side decreases (increases) the 
number density due to annihilation into (production from) other particles, which are 
assumed to be in complete thermal equilibrium. 

It is useful to express the above Boltzmann equation in terms of the dimensionless 
quantities Y-^ = n^/ s and Y^^eq = n^,eq/s. The entropy density is given by s = 
(27r^/45)(7*T^, with g^, being the effective number of relativistic degrees of freedom 
and T the temperature of the universe. We also introduce the dimensionless ratio of 
X mass to the temperature, x = m^/T. Assuming an adiabatic expansion of the 
universe, Eq.(l2]) can then be rewritten as [2j 



The generic picture of x decoupling from the thermal bath is as follows. After 
inflation, the universe becomes radiation-dominated with expansion rate 



(2) 



dt 




(3) 



where Mpi = 2.4 x 10^^ GeV is the reduced Planck mass. The reheat temperate is 
assumed to be high enough for x particles to reach full thermal (chemical as well as 
kinetic) equilibrium]*] Thermal equilibrium is maintained as long as the interaction 
rate F = n^{av) is larger than the Hubble expansion rate H. As the temperature 
decreases, the interaction rate decreases more rapidly than the expansion rate does. 
When the interaction rate falls below the expansion rate, x is no longer kept in ther- 
mal equilibrium and the comoving number density becomes essentially constant. 
This transition temperature is referred to as the freeze-out temperature Tp. 

Analytical expressions for the resulting x relic density are known for the cases 
where decoupling happens when x is non-relativistic {xp = m^/Tp ^ 3) or rela- 
tivistic {xp <^ 3). We discuss these two limiting cases in the following Subsections. 



2.1 Relativist ic case 

First, consider the case where particles x decouple when they are ultra-relativistic 
{xp <^ 3). In this case the equilibrium number density to entropy ratio Y^^cq{x) 
depends on the temperature only through the number of degrees of freedom of 
the thermal bath. Therefore, the final relic abundance is to very good approximation 
equal to its equilibrium value at the time of decoupling: 

>^x,oo = Yxi^ ^ oo) = K^,eq(a;F) = 0.28 iges/9*ixp)) , (5) 

where 

{Oy (for bosons) , , , 

3(7^/4 (for fermions) , 

with being the number of internal (e.g., spin or color) degrees of freedom of 
X- Following the conventional notation, we express the x '^elic density as Q^. = 
m^SoY^^oo/ Pc, where pc = SH^Mp^ is the present critical density of the universe, and 
So — 2900 cm~^ is the present entropy density. This yields 

,2 _ ^7 o in-2 9cs ( rrix 



^ g,{xp) Vl eV/ 

It should be noted that the relic density is simply proportional to the mass of the 
particle in the relativistic case. 



*The case where the reheat temperature is too low for x to attain chemical equilibrium has been 
discussed in [TTl [121 . 
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2.2 Non— relativist ic case 



For the non-rclativistic case where xp ^ 3, the rehc abundance strongly depends 
on the freeze-out temperature Tp because the equihbrium abundance Yy,^cqix) is ex- 
ponentially suppressed as the temperature decreases. The temperature dependence 
of the thermal average of the annihilation cross section is obtained using the Taylor 
expansion in powers of the velocity squared: 

(av) = a + b{v^) + 0{{v^)) = a + - + ( ^] . (8) 



X \x'^ 



The numerically-evaluated correct relic abundance is reproduced with an accuracy 
of a few percent using the approximate analytic formula 

K^,oo = Yxi^ ^ oo) = ^ — . (9) 

1.3 m^Mpl^y g^:(xF){a/xF + 3b/ Xp) 

For WIMPs with electroweak scale mass, freeze-out occurs at ~ 20. The corre- 
sponding scaled relic density is then given by 

Note that the relic density of a non-relativistic particle is inversely proportional to 
its annihilation cross section, but does not depend explicitly on its mass. 



3 Abundance of semi— relativist ically decoupling 
particles 

In the previous Section, we reviewed the known relativistic and non-relativistic ap- 
proximate formulae for the relic abundance. The main aim of this Section is to find 
a simple method applicable between the two regimes: an analytic estimate of the 
relic abundance of semi-relativistically decoupling particles {xp 3). 

One of the key quantities that determine the freeze-out temperature is the 
thermally-averaged cross section. In the non-relativistic case, expressions for the 
equilibrium number density as well as for the thermally-averaged cross section times 
velocity are rather simple. For a semi-relativistic particle, however, the thermally- 
averaged cross section involves multiple integrals and cannot be expanded with re- 
spect to the velocity nor to the mass. Here we discuss a method of approximating 
the thermally-averaged cross section, by interpolating between its relativistically 
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and non-relativistically expanded expressions. We employ the Maxwell-Boltzmann 
distribution for the equihbrium number density [7]: 

The thermal average of the cross section is then obtained as [7j 

1 f°° 

where Ki{x) and K2{x) is the first and second modified Bessel function of the second 
kind; some properties of these functions are given in Appendix A. 

At first glance the use of the Maxwell-Boltzmann distribution seems improper 
for particles that are semi-relativistic at decoupling, let alone for ultra-relativistic 
particles. However, we will argue in Appendix B that this should still yield accurate 
results for the final relic density. This is partly due to cancellations between the 
numerator and denominator of Eq. (fT2l) . and partly due to the fact that the final 
result for f2v becomes less sensitive to decreases. 

As examples of the annihilation of particles, we consider the pair annihilation 
processes of Dirac and Majorana fermions (e.g. neutrinos) into a pair of massless 
fermions, XX ff [H [131 E] • It should, however, be noticed that this assumption 
includes more general cases of any other species annihilating from 5*— or P— wave 
initial states. In a renormalizable model, the annihilation is mediated by some heavy 
particle: for example, a Z-boson with tiny coupling with x, or a new spin-1 boson U 
|15j . We assume that the mass of this exchange particle is much larger than rriy,, so 
that the annihilation amplitude can be described through an effective four-fermion 
interaction. 

The annihilation of two Dirac fermions proceeds from an S*— wave, and the re- 
sulting cross section can be parameterized as 

^sv = Y^, (13) 

where s is the center-of-mass energy squared. G denotes the coupling constant of 
the four-fermion interaction (e.g. the Fermi coupling constant, Gp = 1.17 x 10~^ 
GeV^^). Finally, v is the relative velocity defined as 



(pa ■ PbY - 
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where Pa,b and Ea,b are the four-momenta and energies of the two incident particles 
labeled A and B. The resulting thermally averaged cross section is given by 

y2 



/■OO 

/ dte{e + x'fKi{2V^^^T^). (15) 



Its relativistic and non-relativistic limits read 

respectively. A general expression for {asv) should reproduce these results for x — > 
and x — s> oo, respectively. A simple possibility is 

C^ml (VI 5 + 4x 



It should be noticed that this choice is not unique. 

Let us turn to the case of the annihilation from a P— wave initial state, which is 
e.g. true if x is a Majorana fermion. Eq. flT^ should then be replaced by 

Thermal averaging leads to 

{(Jpv) 



2567rm4TK| 



/•oo 

/ ds sis-AmlfKii^/Ts/T) 

POO 

/ dt t\f + x^)Ki{2VWT^) . (19) 
Jo 



Following the same steps as in the 5— wave case, we find the following interpolation: 

G'^ml /12 3 + 6x 
167r \x'^ [1 + x] 



apv).,, = -^[-,+ J—-T, . (20) 



Figure [Ushows the ratio of the approximate to the exact cross section {crv).,,/ {crv) 
for the S— (solid line) and P— wave (dashed) cases. Note that this ratio depends 
only on x. We see that our approximate expressions reproduce the exact ones with 
accuracy of better than 2% (0.5%) for annihilation from an S— (P— )wave, even in 
the semi-relativistic region (x ~ 3). 

Using Eqs. ffT7|) or fl^Ul) instead of the exact expressions ffT^ or f[T^ greatly reduces 
the numerical effort required to solve the Boltzmann equation ([3]). At the cost of 
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Figure 1: Ratio of the approximate to the exact thermally averaged annihilation cross sections 
(cv)app/(o'v) as a function of x for annihilation from an S— (solid line) and P— wave (dashed) 
initial state. 

some further loss of accuracy, an even faster estimate of the rehc density can be 
obtained by using an approximate solution of the Boltzmann equation instead of 
the accurate evaluation of the relic abundance by solving the Boltzmann equation 
numerically. The determination of the temperature where some interaction decouples 
plays an important role in the analytical prediction of the relic abundance. Indeed, 
for non-relativistically decoupling particles the relic abundance is sensitive to the 
freeze-out temperature xp because of the Boltzmann suppression. In this case, a 
rough estimate of Xi;' obtained by equalizing the interaction rate F and the Hubble 
expansion rate H is not sufficient to make an accurate prediction. Here we show that 
for semi-relativistically decoupling particles a simple comparison of the interaction 
rate and the Hubble expansion rate still gives a reasonably accurate result for the 
final X relic density. 

We define the freeze-out temperature by equalizing the interaction rate and the 
Hubble expansion ratejl 



We then simply assume that does not change after decoupling from the thermal 

''^This definition of the freeze-out temperature should not be used for the prediction of the relic 
abundance through Eqs.© and (|10p . 




(21) 
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Figure 2: Scaled freeze-out temperature xp (right) and scaled relic abundance fl^^h^ (left) as 
function of TO;^. In the right (left) frame the lower (upper) curves are for Majorana fermions, 
and the upper (lower) curves for Dirac fermions. In the left frame the solid curves show exact 
solutions of Eq.(l3]), and the dashed curves our analytic approximations. The dotted curve shows 
the non-relativistic approximation, Eq. (|10p . for P— wave annihilation. Here we take G = Gp = 
1.17 X 10-5 GeV^ = 10 and = 2. 



bath, so that 



Let us see to what extent this method can reproduce the correct rehc abundance. 
As an example, we consider the pair annihilation of neutrino-like particles via the 
mediation of the weak SM gauge bosonj^]. In the left frame of Fig. [2] we plot the 
relic abundance f^^^^ function of m^. The solid curves show predictions for the 
relic abundance obtained by solving the Boltzmann equation ([3]) numerically, while 
the dashed curves have been obtained using the analytic approximation described 
above. The upper curves are for Majorana fermions annihilating from a P— wave, 
and the lower curves are for Dirac fermions case annihilating from an 5— wave. We 



take G = Gf = 1.17 x 10"^ GeV^ g, = 10 and = J. The right frame shows the 
corresponding values of xp. 



^Our primary concern here is to test our approximation for the relic abundance, and thus, as 
an illustration, we take such an unrealistic setup. 

^ Since we concentrate on the ratio of the exact to the approximate relic densities, we discard 
the temperature dependence of g*, which is basically an overall factor. 
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This figure shows that our very simple analytical treatment reproduces the correct 
relic density with an error of at most 20% (5%) for semi-relativistically decoupling 
particles annihilating from an S— (P— )wave. Not surprisingly, our treatment be- 
comes exact for particles that are relativistic at decoupling^ In the Dirac case, our 
approximation coincides with the exact relic abundance for = 1 GeV, corre- 
sponding to a;i;' ~ 16; however, the deviation becomes larger again for larger m^. 
Therefore, our method is not applicable for the entire region of cold relics. Instead, 
one could switch to the usual non-relativistic treatment described in Sec. 2.2 at the 
cross-over value, i.e. at xp = 16. In the P— wave case the cross-over already occurs 
at rri)^ = 30 MeV, corresponding to Xp ~ 4.5. The dotted curve shows that the 
non-relativistic approximation is already quite reliable at this point. We can thus 
smoothly match our approximation to the usual non-relativistic treatment for both 
S— and P— wave annihilation. 



4 Semi— relativistic dark matter? 

As a first application, let us analyze whether semi-relativistically decoupled particles 
(xi? ~ 3) can be a dark matter candidate, whose cosmological abundance should be 
fi^/i^ < 0.1. The final number density of such a particle is of order of Yy.^eq{x — 3) ~ 
10^^. Combining this value with the observed amount of dark matter, the upper 
bound of the semi-relativistic particle turns out to be < 100 eV, which would 
thus decouple at temperatures of a few dozen eV. These particles would therefore still 
be ultra-relativistic during the formation of ^He. Moreover, the effective coupling 
would have to be very large, G ~ 10'^ GeV~^. Such a scenario is therefore tightly 
constrained. 

Note that a x particle with ~ 100 eV could only annihilate into light neutri- 
noso Moreover, the exchange particle would also need to be quite light, with mass 
< 30 MeV if all couplings are < 1. 

The simplest case is that of a real scalar x- Since it only adds a single degree of 
freedom, its presence would not be in serious conflict with current BBN constraints 



^We will see in Appendix B that one should use the correct Fermi-Dirac or Bose-Einstein 

expression for n^Mq in Eq. ()22p in order to accurately predict the relic density for xp ^ 1- 

*The neutrinos themselves aren't in thermal equilibrium with the photons any more at T < 100 
eV. However, as long as T S> rrii, they would still have a thermal distribution. The Boltzmann 
equation ((2]) should therefore still be applicable, if T is taken to be the temperature of the neutrinos, 
which is somewhat lower than the photon temperature. 
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|16j . In a renormalizable theory, XX ~^ could then proceed either through ex- 
change of a fermion in the t— or m— channel, or through boson exchange in the 
s— channel. 

In the former case, the exchange particle would have to be an SU{2) doublet, 
if the low-energy theory only contains left-handed, SU (2) doublet neutrinos. The 
presence of such a light SU (2) doublet fermion is excluded by LEP data. In principle 
the light neutrinos might also be Dirac particles, allowing xX ~^ ^bPr annihilation 
via exchange of a singlet fermion, possibly u^i itself. However, one would then need 
uji to also be in thermal equilibrium, increasing the number of additional degrees of 
freedom present during BBN to an unacceptable level. 

For a real scalar x, s— channel exchange could only proceed through another 
scalar 0. However, a scalar can only couple to z/^zT^ or its hermitean conjugate. 
This scenario would therefore again require uji to have been in equilibrium. We 
therefore conclude that such a light x particle cannot be a real scalar!^ 

If X is a complex field, one needs = 2, which is only marginally compatible 
with BBN [16]. In principle, xX ~^ could then proceed through t— or m— channel 
exchange of an SU{2) singlet. However, then either x or this light exchange particle 
would have to carry hypercharge, so that it would have been produced copiously in 
Z decays. The argument against s— channel exchange of a scalar is the same as for 
real x- 

However, a complex x, either a scalar or a Weyl fermioiJi, could couple to a 
new gauge boson U, which in turn could couple to z^lI^. This new boson would 
contribute three additional bosonic degrees of freedom at T > niu. Since we already 
added two degrees of freedom in x, consistency with BBN would require mu ^ 1 
MeV. This in turn would require the coupling of U to left-handed neutrinos to exceed 
0.01, assuming its coupling to x is perturbative. By SU (2) invariance, U would have 
to couple with equal strength to left-handed charged lepton. This combination of 
U boson mass and coupling is excluded, by a large margin, for both the electron 
and muon family by measurements of the respective magnetic moments [15]. No 
analogous measurement exists for the third generation, so a t/— boson with few MeV 
mass coupling exclusively to third generation leptons and x particles might still be 



^ Of course, this argument does not exclude the possibihty that a much heavier real singlet scalar 

X could be cold Dark Matter [17] . 

massive two-component Weyl fermion can equivalently be described by a four-component 

Majorana fermion. 
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compatible with laboratory data. However, given that /i and r neutrinos are known 
to mix strongly [18], a gauge invariant model where U couples to but does not 
couple to muons is difficult, if not impossible, to construct. 

Finally, such a light x particle would form hot, or at least warm. Dark Matter. 
This possibility is strongly constrained by observations of early structures in the 
universe, in particular the "Lyman— a forest" [19j. Such a x particle could thus at 
best form a sub-dominant component of the total Dark Matter. 

In combination, these arguments strongly indicate that semi-relativistically de- 
coupling particles should not be absolutely stable. In the next Section we will show 
that such particles may nevertheless have a role to play in the history of the Universe, 
if they are metastable. 

5 Entropy production by decaying particles 

In this section we demonstrate that semi-relativistically decoupling particles can be 
useful for producing a large amount of entropy, which could dilute the density of 
other relics to an acceptable level. Examples of such relics are decaying gravitinos, 
which can lead to problems with Big Bang nucleosynthesis [20], or supersymmetric 
neutralinos, whose relic density often exceeds the required Dark Matter density by 
one or two orders of magnitude [2T]. The density of such relics will be diluted 
only if the entropy is released after they decouple from the thermal bath. This will 
simultaneously dilute any pre-existing baryon asymmetry. One thus either has to 
increase the efficiency of early baryogenesis, or introduce late baryogenesis after the 
release of the additional entropy. Both possibilities can be realized in the framework 
of Affieck-Dine baryogenesis |22j . 

Generally ^ \9\, out-of-equilibrium decays of long-lived particles can only pro- 
duce a significant amount of entropy if the decaying particle dominates the energy 
density of the Universe prior to its decay. The abundance of non-relativistically 
decoupling particles is suppressed by a factor e~^^, hence their contribution to the 
energy density is small at decoupling. However, after decoupling their energy den- 
sity only drops like oc T^, while that of the dominant radiation component 
decreases like as the Universe cools offi Thermally produced particles can there- 
fore dominate the energy density of the Universe only at temperature T <C e~^^Tp. 
Significant entropy production by the late decay of nonrelativistically decoupling 
particles is therefore only possible if they are simultaneously very massive and quite 
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long-lived. For semi-relativistic particles, on the other hand, the abundance at de- 
coupling is large and thus a significant amount of entropy can be produced even 
if their mass is small, since their density will become dominant quite soon after 
decoupling. 

Let us consider the out-of-equilibrium decay of long-lived particles which semi- 
relativistically decoupled from the thermal background. For simplicity we work in 
the instantaneous decay approximation, i.e. we assume that all x particles decay 
at time td = t^, where is the lifetime of x- While this approximation does not 
describe the time dependence of the entropy (or temperature) for t ^ very well, 
it does reproduce the entropy enhancement factor, i.e. the entropy at t ^ r^, quite 
accurately. We assume that x particles were in full thermal equilibrium for suffi- 
ciently high temperatures in the RD epoch. When the temperature decreased to 
T = Tp ^ m^, the x number density froze out. At decoupling, x particles 
contributed a few percent to the total energy density of the universe; however, as 
noted earlier, the ratio of the radiation and x energy densities decreased by a factor 
Td/Tp = a/^f/tx between decoupling and decay of x; here Td refers to the temper- 
ature at time t = r^, just prior to x decay. If tp, the x energy density at the 
time of the x decay is well approximated by py^ d = ^x^x4^ dominated over the 
radiation. In this case, the ratio of the final entropy density Sf after the decay to 
the initial entropy density Si before the decay is given by [8] 

^ = 0.82^y^^^^^i^, (23) 

for Si <^ Sf. Here Y^d = n^,d/ Si is proportional to the x abundance just prior to its 
decay. 

In the light of the BBN prediction of the primordial abundances of the light 
elements, the x lifetime is constrained as < 1 sec [23] • Equations ffTUl) and fl2^ 
show that the entropy ratio is proportional to the relic density ^^^"^ that x would 
have if it were stable. We saw in Fig. 2 that for fixed coupling G this quantity 
is maximal if Tp ~ m^; more accurately, the maximum of ^^^"^ achieved for 
xp = 1.8 (2.1) if X particles annihilate from an S— (P— )wave initial state. Entropy 
production by late x decays is thus most efficient when the x particles decoupled 
semi-relativistically, with their lifetime fixed to the maximal value of ~ 1 sec. 

We can construct a feasible scenario that fulfills these conditions by introduc- 
ing a sterile neutrino which mixes with an ordinary neutrino. Here we treat both 
m-^ and the mixing angle 6 as free parameters. In sharp contrast to conventional 
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cosmological scenarios with sterile neutrinos [lO], the sterile neutrino is assumed to 
be in thermal equilibrium in the early universe. In ordinary sterile neutrino mod- 
els, thermal equilibrium is not reached because the Yukawa coupling of the sterile 
neutrino with SM particles is tiny. One possible method for the x Pci^ir production 
and annihilation to reach thermal equilibrium is to extend sterile neutrino models 
by adding another hypothetical boson Z' . Let it have coupling gz'f with the SM 
fermion pair //, and gz'x with the sterile neutrino pair. If the Z'-boson mass m^/ 
is larger than the x energy, the x annihilation cross section has the form of Eq. ( fTSl) 
with G = gz'x9z'f/fn'z'- Although gz'f and rriz' are constrained by high energy ex- 
periments, gz'x can be as large as unity. Therefore, x annihilation can be in thermal 
equilibrium before its semi-relativistic decoupling. Decoupling occurred at T ~ 
if ~ 1 GeV ■ (3 ■ lO'^ GeV-ycY^^gl^^. 

In order to estimate the amount of entropy released by the decay of sterile neu- 
trinos in this setup, we have to calculate their lifetime. For simplicity we ignore 
propagator effects. When the sterile neutrino mass is smaller than the IV— boson 
mass mw = 80 GeV, it decays into three SM fermions, with decay width 



27 - 16 sin^ Ow + — sin^ Ow 
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(24) 



where sin^ 6w = 0.23 is the weak mixing angle. When the sterile neutrino mass is 
larger than the Z— boson mass mz = 91 GeV, the sterile neutrino predominantly 
decays into a SM gauge boson and a lepton. Its decay width is then proportional to 



m: 



X' 



and given by 
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w 



1 + 



2m|. 



^ 1 + 



In the in-between case where m\Y < < mz, we obtain 



2m| 



sin^^. (25) 



+ 



2^^ 

8V27T 

11 



m 



w 



1 + 



2m^ 



20 sin^ ew + — sin^ Ow 



80 

y 



1927r3 



sin^ e . 



(26) 



Figure [3] shows contours of the entropy increase Sf/si due to sterile neutrino 
decay in the (1/x/G, sin^) plane. We set the freeze-out temperature to xp = 2.1, 
which maximizes m^Y^ i] this can be achieved by chosing the mass appropriately. 
The thick line indicates the BBN limit on the sterile neutrino lifetime, = 1 sec. 
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Figure 3: Contour of the entropy increase s//si caused by semi-relativistic sterile neutrino decay 
in the {l/VG,sm9) plane. We choose such that xp = 2.1. The solid line indicates the BBN 
limit on the sterile neutrino lifetime = 1 sec. 

Eg. (1231) shows that for given neutrino mass, the released entropy will be maximal if 
6 is chosen such that r reaches this upper limit. 

The behavior of the contours in Fig. [3] is easy to understand from Eg. (1231) . In 
the relevant limit ^ 1 and keeping (7* constant, we have oc 6~'^m~^ {6~'^m~^) 
for < (>) mw- The entropy ratio thus scales as O^^niy^^^'^ oc 9^^G {O^^rriy^^^'^ oc 
9^^G^^^) for < (>) mw- Along the = 1 sec contour, the entropy release 
increases proportional to m^^ oc G^^/^ both for < m^y and for m^^ > mw Fig. [3] 
can be extended to even smaller G, i.e. larger Z' masses, so long as m^ is smaller 
than the re-heat temperature after inflation, so that x "was in thermal eguilibrium 
in the RD epoch. If at the same time 6 is decreased so that = 1 sec remains 
constant, very large entropy dilution factors could be realized, 



This result is only valid if the mixing-induced interactions of x are not in thermal 
eguilibrium for T < m^. Since these interactions are also responsible for x decay, 
this assumption is satisfied whenever tp; we saw in the discussion of Eg.( !23|) 




(27) 
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that this strong inequahty is in any case a condition for significant entropy release 
from X decay. 

We finally note that for given the entropy released in x decays is maximal 
if G is so small that x was ultra-relativistic at decoupling, since this maximizes 
Y^^eqixp)- Again setting = 1 sec by appropriate choice of 9, this yields 

££ < 104 . /28) 
Si - 103 GeV ^ ^ 

6 Conclusion 

In this paper we have developed an approximate analytic method for calculating 
the thermally-averaged annihilation cross section of semi-relativistically decoupling 
particles and for estimating their relic density. We have shown that this approximate 
solution can be smoothly matched to the well-known non-relativistic approximation 
at the point of intersection. We have argued that such relics cannot form the observed 
cosmological dark matter. However, we pointed out that the late decay of metastable 
semi-relativistically decoupling relics can be an efficient source of entropy produc- 
tion. As an example of this entropy production mechanism we discussed a scenario 
with a sterile neutrino, and illustrated to what extent entropy can be increased. 
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Appendix A: Modified Bessel Functions 

In this appendix, we summarize some properties of the modified Bessel function. 
Using an integral representation, the modified Bessel function of the second kind is 
defined by 

K.{z) = r'(!+V2) ly " ^'^^^ ^ ^'^^^ ^ ° ■ ^^^^ 
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In particular, the calculation of the relic abundance involves Ki{z) and K2{z), 

/oo 
dt e-"*(t2 - 1)1/2^ Re(z)>0, 

2 roo 

K2{z) = — dt e-^*(t2 - 1)3/2, Ke(z)>0. (30) 
3 Jl 

The lower order terms of the series expansion of Ki{z) and K2{z) are given by 

Ki{z) = - + ■■■, 
z 

K2{z) = ^2-\ + ---- (31) 
The asymptotic expansion of Ky{z) is given by 

/CW~/|e-(l + i^ + ..■). (32) 

Appendix B: Validity of the Maxwell— Bolt zmann 
Distribution 

In the calculations of this paper we used the Maxwell-Boltzmann (MB) distribution 
also for particles that were semi-relativistic at decoupling; this assumption is e.g. 
implicit in Eq. f|T2l) . At first sight this seems quite dangerous. For example, at T = 
m^, i.e. a; = 1, the MB result for ?T-x,eq overestimates the Fermi-Dirac distribution by 
about 7%, and underestimates the Bose-Einstein distribution by about 10%. Since 
X annihilation always involves two x particles, one might assume that the total 
error associated with the use of the MB distribution is about twice as large. In this 
Appendix we show that the MB distribution can indeed be used to compute the 
thermally averaged cross section and the decoupling temperature as long as x^? > 1. 
For smaller xp-, one has to use the proper Fermi-Dirac or Bose-Einstein distribution 
only in the very last step, when calculating yj^,oo- 

We begin by expanding the true distribution function, 

kAEd = ^^^7^ - e-x/- (1 ^ e-^^/-) , (33) 

where the upper (lower) sign is for fermionic (bosonic) x particles. Note that the 
correction term in parentheses has exactly the same form as the "statistics factors" 
appearing in the collision term of the full Boltzmann equation Pj. For consistency 
these statistics factors therefore also have to be included. Up to first order in these 
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correction factors, the temperature dependent terms in the integrand defining the 
colhsion term for XX ^ ff processes then read for fermionic /: 

_ (1 _ e-Ef/T _ ^-EjIT ^ c^^-E^^/T ^ c^e-^«/^)] 
= ^-{E^,+E^,)/T . ^(^2 _ 1) ^1 ^ ^-E^JT ^ ^-E^JT _ ^-E,/T _ ^-E^/T^^ 

±(c^_l)(e-^xi/T + e-i?x./T)] . (34) 

here = /x//x,eq is independent of energy as long as x is in kinetic equihbrium 
(through elastic scattering on SM particles); in that case we can equivalently write 

= n^/n^^eq- In order to derive the full collision term, I has to be multiplied with 
the squared matrix element and integrated over phase space [2]. 

In the usual treatment of WIMP decoupling, all the exponential terms in the 
square parentheses are neglected, so that the collision term becomes proportional to 

~ '^x.eq ti^6s the thermally averaged cross section defined in Eq. (fT2|) . Unfortu- 
nately the full correction term introduces additional dependence on the final state 
energies Ef and Ej. In order to keep the numerics manageable, we assume that 
they can be replaced by E^-^ and E^^, respectively. This is certainly true (by energy 
conservation) for the sum Ef + Ej; this has already been used in deriving Eq. (1341) . 
Note furthermore that we'll need the collision term for temperatures > Tp, where 
|cj(, — 1| ^ 1, so that we can approximate — 1 ~ (c^ — l)/2. These approximations 
yield 

J ^ (4 - 1) e-(^«+^«)/^ (e-^-i/^ + e-^«/^)] , (35) 

where k = 1/2 (3/2) for bosonic (fermionic) x particles. In the following we will 
assume x to be fermionic, which according to Eq. (135|) should lead to larger deviations 
from the MB result. 

Inserting this corrected collision term into the Boltzmann equation, and following 
the formalism of [7], finally yields a modified thermally averaged cross section times 
initial state velocity: 

-.2 



{av) - 



<eq8(27r)' 



j dE+dE^ds{a ■ F)(s)e-^+/^ 

. [1 _ ^ (g-(i^++i?-)/(2T) ^ ^-iE^-E.)/i2T)^^ ^ (3g) 

with F = 2s^l — 4m^/s, E+ = E^^-^ + £"^3 &nd = E^-^ — E^^. This reduces to 
Eq.f ll2l) if the expression in square parentheses is simply replaced by 1. In the fol- 
lowing we assume that x particles annihilate from an S"— wave. P— wave annihilation 
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would favor larger energies, where the correction terms in Eq. ( l36i) are smaller. Note 
that we also have to use the expanded form (!33l) of the distribution function when 
calculating n^ eq in Eq. (l36|) : otherwise the solution of the Boltzmann equation will 
not yield ~ ^x.eq; including the correction terms, aX. T ^ Tp. 




Figure 4: Various approximations for the thermally averaged cross section as function of the 
scaled inverse temperature x — rn^/T for fcrmionic particles annihilating from an S"— wave. The 
solid (black) curve shows the ratio of the corrected cross section ([36|) to the Maxwell-Boltzmann 
(MB) result (fT2|) . while the dashed (red) curve shows this ratio if Eq. (fl2|) is replaced by our 
approximation pT|) . The dotted (blue) curve is the same as the solid curve in Fig. 1. 



The size of the correction terms in Eq. (l36|) is illustrated by the solid (black) curve 
in Fig.m We see that the correction amounts to less than 2% for all x > 1. This is due 
to a strong cancellation between the corrections in the integrand of Eq. ( !36l) and those 
in the overall factor 1/?t,^„„. The dashed (red) curve shows that for x ~ 2 the errors 
due to the use of the MB distribution and due to our simple parameterization (|T7|) 
add up, leading to a total error of about 2.7% at most. The Fermi-Dirac corrections 
to the thermally averaged cross section begin to be significant for x < 0.5. However, 
here one enters the ultrarelativistic regime, where the final relic density is no longer 
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sensitive to the decoupling temperature. We therefore expect the effect of using the 
MB distribution in Eq. (fT2!) on the final prediction of the relic density to be quite 
small throughout. 
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Figure 5: Effect of using the Maxwell-Boltzmann distribution on the predicted relic density, 
calculated using the approximation Y^^oc = Y-^^^cqixp), for fermionic x particles annihilating from an 
iS— wave initial state. The dashed (red) curve shows the ratio of the prediction using the Maxwell- 
Boltzmann distribution everywhere to the corrected prediction based on Eq. (|33p and (|36p . The 
solid (black) curve shows the analogous ratio, where correct Fermi-Dirac distribution has been 
used to evaluate Y^^eq{xF), but (av) and xp have still been obtained using the MB distribution. 
Parameters are as in Fig. 2. 



This is illustrated in Fig. [5l where the relic density has been calculated using 
the simple assumption Y^^oo = ^x.eq(^i^)j have used the same parameters as in 
Fig. 2. The dashed (red) curve shows that using the MB distribution everywhere will 
overestimate the relic density for < 5 MeV, i.e. for xp ^ I- However, the black 
curve shows that this can easily be corrected by using the Fermi-Dirac distribution 
only in the final step, i.e. when calculating Y^^eq^xp)', {crv) and xp can still been 
calculated using the MB distribution. This validates our treatment in the main text. 
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